/ BBS in a Box 12 / BBS in a box XII-1.iso / Files / Education / M / MathPad 2.15.sit / MathPad 2.15 / examples / curve fit

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 00 09 63 75 72 76 65 20 | 66 69 74 00 00 00 00 00 |..curve |fit.....|
|00000010| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 54 45 58 54 4d 50 41 | 44 00 00 00 00 00 00 00 |.TEXTMPA|D.......|
|00000050| 00 00 00 00 00 04 22 00 | 00 01 56 a8 4b e0 95 a8 |......".|..V.K...|
|00000060| 4b e0 95 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |K.......|........|
|00000070| 00 00 00 00 00 00 00 00 | 00 00 81 81 e5 36 00 00 |........|.....6..|
|00000080| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 67 65 6e 65 72 61 |--------|--genera|
|00000090| 6c 20 6c 65 61 73 74 20 | 73 71 75 61 72 65 73 20 |l least |squares |
|000000a0| 66 69 74 2d 2d 2d 2d 2d | 2d 0d 2d 2d 20 66 69 6e |fit-----|-.-- fin|
|000000b0| 64 73 20 63 6f 65 66 66 | 69 63 69 65 6e 74 73 20 |ds coeff|icients |
|000000c0| 61 5b 69 5d 20 66 6f 72 | 20 62 65 73 74 20 66 69 |a[i] for| best fi|
|000000d0| 74 20 6f 66 20 77 65 69 | 67 68 74 65 64 20 73 75 |t of wei|ghted su|
|000000e0| 6d 20 6f 66 20 66 75 6e | 63 74 69 6f 6e 73 20 6f |m of fun|ctions o|
|000000f0| 66 20 78 0d 0d 66 69 74 | 28 78 29 20 3d 20 73 75 |f x..fit|(x) = su|
|00000100| 6d 28 61 5b 69 2c 31 5d | 2a 66 28 78 29 5b 69 5d |m(a[i,1]|*f(x)[i]|
|00000110| 2c 69 2c 31 2c 6e 70 61 | 72 6d 73 29 0d 0d 66 28 |,i,1,npa|rms)..f(|
|00000120| 78 29 5b 6a 5d 3d 78 5e | 28 6a 2d 31 29 20 20 2d |x)[j]=x^|(j-1)  -|
|00000130| 2d 20 70 6f 6c 79 6e 6f | 6d 69 61 6c 20 61 31 2b |- polyno|mial a1+|
|00000140| 61 32 2a 78 2b 61 33 2a | 78 5e 32 2e 2e 2e 0d 0d |a2*x+a3*|x^2.....|
|00000150| 50 5b 69 2c 6a 5d 3d 66 | 28 78 5b 69 5d 29 5b 6a |P[i,j]=f|(x[i])[j|
|00000160| 5d 20 64 69 6d 5b 6e 64 | 61 74 61 2c 6e 70 61 72 |] dim[nd|ata,npar|
|00000170| 6d 73 5d 0d 50 74 50 3a | 3d 6d 75 6c 74 69 70 6c |ms].PtP:|=multipl|
|00000180| 79 28 74 72 61 6e 73 70 | 6f 73 65 28 50 29 2c 50 |y(transp|ose(P),P|
|00000190| 29 3a 0d 50 74 79 3a 3d | 6d 75 6c 74 69 70 6c 79 |):.Pty:=|multiply|
|000001a0| 28 74 72 61 6e 73 70 6f | 73 65 28 50 29 2c 79 29 |(transpo|se(P),y)|
|000001b0| 3a 0d 69 6e 76 3a 3d 69 | 6e 76 65 72 74 28 50 74 |:.inv:=i|nvert(Pt|
|000001c0| 50 29 3a 0d 61 3a 3d 6d | 75 6c 74 69 70 6c 79 28 |P):.a:=m|ultiply(|
|000001d0| 69 6e 76 2c 50 74 79 29 | 3a 0d 0d 64 61 74 61 3d |inv,Pty)|:..data=|
|000001e0| 72 65 61 64 28 78 79 64 | 61 74 61 29 0d 78 5b 69 |read(xyd|ata).x[i|
|000001f0| 5d 3d 64 61 74 61 5b 69 | 2c 31 5d 3b 20 79 5b 69 |]=data[i|,1]; y[i|
|00000200| 5d 3d 64 61 74 61 5b 69 | 2c 32 5d 20 64 69 6d 5b |]=data[i|,2] dim[|
|00000210| 6e 64 61 74 61 5d 3b 20 | 0d 0d 6e 64 61 74 61 3d |ndata]; |..ndata=|
|00000220| 63 6f 75 6e 74 28 64 61 | 74 61 29 0d 6e 70 61 72 |count(da|ta).npar|
|00000230| 6d 73 3d 33 0d 0d 61 3a | 7b 7b 35 2e 35 38 32 7d |ms=3..a:|{{5.582}|
|00000240| 2c 7b 30 2e 31 35 39 7d | 2c 7b 30 2e 30 30 36 7d |,{0.159}|,{0.006}|
|00000250| 7d 0d 0d 70 6c 6f 74 20 | 64 61 74 61 0d 70 6c 6f |}..plot |data.plo|
|00000260| 74 20 66 69 74 28 58 29 | 0d 0d 2d 2d 2d 2d 2d 2d |t fit(X)|..------|
|00000270| 2d 2d 2d 2d 20 6d 61 74 | 72 69 78 20 6f 70 65 72 |---- mat|rix oper|
|00000280| 61 74 69 6f 6e 73 20 2d | 2d 2d 2d 2d 2d 2d 2d 2d |ations -|--------|
|00000290| 2d 2d 2d 0d 6d 75 6c 74 | 69 70 6c 79 28 41 2c 42 |---.mult|iply(A,B|
|000002a0| 29 5b 69 2c 6a 5d 20 3d | 20 73 75 6d 28 41 5b 69 |)[i,j] =| sum(A[i|
|000002b0| 2c 6b 5d 2a 42 5b 6b 2c | 6a 5d 2c 6b 2c 0d 20 20 |,k]*B[k,|j],k,.  |
|000002c0| 20 20 20 20 20 20 20 31 | 2c 63 6f 75 6e 74 28 42 |       1|,count(B|
|000002d0| 29 29 20 64 69 6d 5b 63 | 6f 75 6e 74 28 41 29 2c |)) dim[c|ount(A),|
|000002e0| 63 6f 75 6e 74 28 42 5b | 31 5d 29 5d 0d 74 72 61 |count(B[|1])].tra|
|000002f0| 6e 73 70 6f 73 65 28 41 | 29 5b 69 2c 6a 5d 20 3d |nspose(A|)[i,j] =|
|00000300| 20 41 5b 6a 2c 69 5d 20 | 64 69 6d 0d 20 20 20 20 | A[j,i] |dim.    |
|00000310| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|00000320| 20 20 20 20 5b 63 6f 75 | 6e 74 28 41 5b 31 5d 29 |    [cou|nt(A[1])|
|00000330| 2c 63 6f 75 6e 74 28 41 | 29 5d 0d 69 6e 76 65 72 |,count(A|)].inver|
|00000340| 74 28 41 29 20 3d 20 61 | 64 6a 6f 69 6e 74 28 41 |t(A) = a|djoint(A|
|00000350| 29 2f 64 65 74 28 41 29 | 0d 61 64 6a 6f 69 6e 74 |)/det(A)|.adjoint|
|00000360| 28 41 29 20 3d 20 74 72 | 61 6e 73 70 6f 73 65 28 |(A) = tr|anspose(|
|00000370| 63 6f 66 61 63 74 6f 72 | 28 41 29 29 0d 63 6f 66 |cofactor|(A)).cof|
|00000380| 61 63 74 6f 72 28 41 29 | 5b 69 2c 6a 5d 20 3d 20 |actor(A)|[i,j] = |
|00000390| 28 2d 31 29 5e 28 69 2b | 6a 29 2a 0d 20 20 20 20 |(-1)^(i+|j)*.    |
|000003a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 64 |        |       d|
|000003b0| 65 74 28 73 75 62 6d 61 | 74 72 69 78 28 41 2c 69 |et(subma|trix(A,i|
|000003c0| 2c 6a 29 29 0d 73 75 62 | 6d 61 74 72 69 78 28 41 |,j)).sub|matrix(A|
|000003d0| 2c 6b 2c 6c 29 5b 69 2c | 6a 5d 20 3d 20 41 5b 69 |,k,l)[i,|j] = A[i|
|000003e0| 2c 6a 5d 20 77 68 65 6e | 20 69 3c 6b 20 61 6e 64 |,j] when| i<k and|
|000003f0| 20 6a 3c 6c 2c 0d 20 20 | 20 20 20 20 20 20 20 20 | j<l,.  |        |
|00000400| 20 20 41 5b 69 2b 31 2c | 6a 5d 20 77 68 65 6e 20 |  A[i+1,|j] when |
|00000410| 69 b3 6b 20 61 6e 64 20 | 6a 3c 6c 2c 0d 20 20 20 |i.k and |j<l,.   |
|00000420| 20 20 20 20 20 20 20 20 | 20 41 5b 69 2b 31 2c 6a |        | A[i+1,j|
|00000430| 2b 31 5d 20 77 68 65 6e | 20 69 b3 6b 20 61 6e 64 |+1] when| i.k and|
|00000440| 20 6a b3 6c 2c 0d 20 20 | 20 20 20 20 20 20 20 20 | j.l,.  |        |
|00000450| 20 20 41 5b 69 2c 20 20 | 6a 2b 31 5d 20 77 68 65 |  A[i,  |j+1] whe|
|00000460| 6e 20 69 3c 6b 20 61 6e | 64 20 6a b3 6c 20 64 69 |n i<k an|d j.l di|
|00000470| 6d 0d 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |m.      |        |
|00000480| 20 20 20 20 20 20 20 20 | 20 20 5b 63 6f 75 6e 74 |        |  [count|
|00000490| 28 41 29 2d 31 2c 63 6f | 75 6e 74 28 41 29 2d 31 |(A)-1,co|unt(A)-1|
|000004a0| 5d 0d 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |].......|........|
|000004b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000004c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000004d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000004e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000004f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000500| 00 00 01 00 00 00 01 24 | 00 00 00 24 00 00 00 32 |.......$|...$...2|
|00000510| 69 64 65 57 69 6e 64 6f | 77 28 64 69 61 6c 6f 67 |ideWindo|w(dialog|
|00000520| 29 3b 0d 20 20 20 20 20 | 20 6c 6f 6f 70 3d 46 41 |);.     | loop=FA|
|00000530| 09 63 75 72 76 65 20 66 | 69 74 02 00 00 00 50 61 |.curve f|it....Pa|
|00000540| 72 74 53 49 54 78 00 00 | 00 00 00 00 00 00 00 01 |rtSITx..|........|
|00000550| 00 00 50 61 72 74 53 49 | 54 78 00 00 00 00 00 00 |..PartSI|Tx......|
|00000560| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000570| 00 00 a9 7e bd ee 00 00 | 00 00 00 00 01 56 30 4c |...~....|.....V0L|
|00000580| 29 3b 0d 20 20 20 20 20 | 20 0d 20 20 20 20 20 64 |);.     | .     d|
|00000590| 65 66 61 75 6c 74 3a 09 | 2f 2a 20 65 64 69 74 61 |efault:.|/* edita|
|000005a0| 62 6c 65 20 74 65 78 74 | 20 66 69 65 6c 64 73 20 |ble text| fields |
|000005b0| 2a 2f 0d 20 20 20 20 20 | 20 62 72 65 61 6b 3b 0d |*/.     | break;.|
|000005c0| 20 20 20 20 7d 0d 20 20 | 20 7d 0d 20 20 20 72 65 |    }.  | }.   re|
|000005d0| 74 75 72 6e 28 64 69 61 | 6c 6f 67 29 3b 0d 7d 0d |turn(dia|log);.}.|
|000005e0| 0d e6 0d f4 0d f6 0e 20 | 0e 4a 0f 2c 8e 52 0e 52 |....... |.J.,.R.R|
|000005f0| 0e 60 8e 8c 0e 8c 8e 9c | 0e 9c 0e c6 0e e2 8f 08 |.`......|........|
|00000600| 00 00 00 20 00 00 00 03 | 00 02 3f f9 8e fa 35 12 |... ....|..?...5.|
|00000610| 94 e9 c8 ae 01 36 01 25 | 00 05 00 30 00 bd 01 25 |.....6.%|...0...%|
|00000620| 01 3f 00 2f 00 00 01 00 | 00 00 01 24 00 00 00 24 |.?./....|...$...$|
|00000630| 00 00 00 32 00 2b 76 84 | 05 e2 00 00 00 1c 00 32 |...2.+v.|.......2|
|00000640| 00 00 50 52 65 66 00 00 | 00 0a 00 80 ff ff 00 00 |..PRef..|........|
|00000650| 00 00 00 2c 7b 0c 00 00 | 00 00 00 00 00 00 00 00 |...,{...|........|
|00000660| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000670| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
+--------+-------------------------+-------------------------+--------+--------+